In silico and gene expression analysis of the acute inflammatory response of gilthead seabream (Sparus aurata) after subcutaneous administration of carrageenin

Inflammation is one of the main causes of loss of homeostasis at both the systemic and molecular levels. The aim of this study was to investigate in silico the conservation of inflammation-related proteins in the gilthead seabream (Sparus aurata L.). Open reading frames of the selected genes were used as input in the STRING database for protein–protein interaction network analysis, comparing them with other teleost protein sequences. Proteins of the large yellow croaker (Larimichthys crocea L.) presented the highest percentages of identity with the gilthead seabream protein sequence. The gene expression profile of these proteins was then studied in gilthead seabream specimens subcutaneously injected with carrageenin (1%) or phosphate-buffered saline (control) by analyzing skin samples from the injected zone 12 and 24 h after injection. Gene expression analysis indicated that the mechanisms necessary to terminate the inflammatory response to carrageenin and recover skin homeostasis were activated between 12 and 24 h after injection (at the tested dose). The gene analysis performed in this study could contribute to the identification of the main mechanisms of acute inflammatory response and validate the use of carrageenin as an inflammation model to elucidate these mechanisms in fish.


Introduction
Inflammation is a temporary local response characterized by heat, redness, swelling, pain, and functional disorders (Nathan 2002). It is triggered by the innate immune system to protect the host against tissue damage, and its aim is to restore physiological functions when homeostatic mechanisms are insufficient (Medzhitov 2008). Inflammation consists of a sequence of events. Depending on diverse factors (such as the causal agent of the inflammation or the anatomical site where it is triggered), inflammatory mediators (including prostaglandins, histamines, cytokines, reactive oxygen species, and enzymes) are released by tissue-resident immune and non-immune cells to increase the blood flow and cause vasodilation (Calder et al. 2013). Capillary permeability is subsequently increased due to the retraction of endothelial cells, allowing access of many large and soluble molecules to the inflammation site (Calder et al. 2013). The third main event is the recruitment of leucocytes (granulocytes, neutrophils, eosinophils Abstract Inflammation is one of the main causes of loss of homeostasis at both the systemic and molecular levels. The aim of this study was to investigate in silico the conservation of inflammation-related proteins in the gilthead seabream (Sparus aurata L.). Open reading frames of the selected genes were used as input in the STRING database for protein-protein interaction network analysis, comparing them with other teleost protein sequences. Proteins of the large yellow croaker (Larimichthys crocea L.) presented the highest percentages of identity with the gilthead seabream protein sequence. The gene expression profile of these proteins was then studied in gilthead seabream specimens subcutaneously injected with carrageenin (1%) or phosphate-buffered saline (control) by analyzing skin samples from the injected zone 12 and 24 h after injection. Gene expression analysis indicated that the mechanisms necessary to terminate the inflammatory response to carrageenin and recover skin homeostasis were activated between 12 and 24 h after injection (at the tested dose). The gene analysis performed in this study could contribute to the identification of the main mechanisms of acute and basophils, and monocytes-macrophages) from the blood capillaries to the surrounding tissues and subsequently to the inflammation site due to the repertoire of chemokines and adhesion molecules released by endothelial cells (Calder et al. 2013). Finally, inhibitory mechanisms are normally activated to terminate the inflammatory process and initiate the repair of damaged tissue (Calder et al. 2013), thus preventing unnecessary damage and restoring homeostasis (Chen and Nuñez 2010;Henson 2005). However, if these mechanisms fail to restore homeostasis, inflammation may persist, resulting in a pathological state of chronic inflammation and even chronic diseases (Nathan and Ding 2010;Serhan et al. 2011).
Although inflammatory processes have been extensively studied in mammals, they have not been sufficiently studied in fish. The intensive fish production in aquaculture, one of the fastestgrowing sectors in the world in the past few decades, is very frequently associated with inflammation caused injuries, especially wounds and ulcers in the skin of farmed fish (Chen and Nuñez 2010). Moreover, fish are highly susceptible to stress due to the intimate contact with their environment, which can compromise their welfare and promote illness (Delamare-Deboutteville et al. 2006). The final consequences of inflammation and the outbreak of diseases in fish are substantial economic losses (Balcázar et al. 2006;Esteban 2012). Therefore, the search and use of non-invasive methods, such as in silico analyses, can be beneficial for the study of molecular processes that have scarcely been studied in fish, such as inflammation Romano and Tatonetti 2019). Exploring the conservation of pro-inflammatory and antiinflammatory proteins and other inflammationrelated molecules among species, as well as their functions, could contribute to the identification of the patterns and signal transduction pathways that regulate this complex process (Hawiger and Zienkiewicz 2019). Since physical interactions between proteins play an important role in this process, study methods are increasingly based on predicting protein-protein interactions purely from protein sequences to identify proteins that have not been previously studied (Makrodimitris et al. 2020). The public online database Search Tool for the Retrieval of Interacting Genes and Proteins (STRING; https:// string-db. org/) provides accurate information on protein associations (interactions) by combining different sources and techniques, such as gene co-expression and text mining, and can assign a protein interaction probability score to each physical interaction and functional association (Szklarczyk et al. 2019). These scores have been used to reconstruct protein-protein interaction networks (Kummer et al. 2018). However, the complexity of the studied organisms and the information extracted from the database may vary to such a degree that a protein interaction may not coincide with the organism of interest or may not offer contextual information, necessitating experimental validation of the obtained results (Björne and Salakoski 2019).
In this study, we used carrageenin to better understand the mechanism of inflammation. Carrageenin is a high-molecular-weight sulfated mucopolysaccharide derived from the cell wall of the red algae commonly called Irish moss or carrageen moss (Chondrus crispus). It has been used for decades as a model of acute inflammation in rats, mice, and guinea pigs and in in vitro inflammation experiments in humans (Bhattacharyya et al. 2008a, b;Bhattacharyyaa et al. 2010;Levy 1969;Morris 2003;Winter et al. 1962). Carrageenin is composed of 1,3α-1,4β-galactans having one (κ-), two (ι-), or three (λ-) sulfates per disaccharide unit and can produce an inflammatory response associated with edema, hyperalgesia, and erythema and an acute response to thermal and mechanical stimuli. When self-associating into helical structures, κ-and ι-carrageenin can form rigid and flexible gels, respectively, in ionic solutions and seems to be related to the immunostimulant activity against bacterial infections in teleost fish (Cheng et al. 2007;Fujiki et al. 1994Fujiki et al. , 1997Necas and Bartosikova 2013). In contrast, λ-carrageenin does not form helices or gels, but it is the main compound responsible for inducing an acute inflammation response when injected in rodents (Fujiki et al. 1997). The few studies investigating the effects of carrageenin on fish have reported that intraperitoneal injection of carrageenin induces acute inflammation characterized by leucocyte infiltration and significant abdominal edema a few hours after administration, peaking several hours later (Huang et al. 2014;Martins et al. 2006;Matushima and Mariano 1996;Ribas et al. 2016). After being injected, carrageenin seems to remain extracellularly for several hours, causing chronic inflammation until macrophages begin to engulf and store it in cytoplasmic vesicles, eventually removing all the extracellular carrageenin by the end of the inflammation process (Timur et al. 1977).
In light of these findings, the aim of this study was to evaluate the involvement of 40 inflammationrelated proteins (cellular markers and proinflammatory and anti-inflammatory proteins) in the inflammatory process evoked by the subcutaneous administration of carrageenin. Carrageenin was selected because it is known to trigger inflammation in mammals (Levy 1969;Winter et al. 1962;Morris 2003), while the gilthead seabream was chosen as a fish species of interest to marine aquaculture (APROMAR 2019). STRING was used to compare the protein sequence of the genes studied in gilthead seabreams with protein sequences in other similar fish species. The functional associations and interactions between the proteins involved in the same process and their level of conservation were also studied. To our knowledge, this is the first study to evaluate the modulation of gene expression after carrageenin administration in fish.

Protein prediction and functional protein association network analysis
Based on previous inflammation studies on mammals and zebrafish (Campos-Sánchez and Esteban 2021; Chen and Nuñez 2010) and taking into account the published gilthead seabream genome, 40 inflammation-related proteins were selected to investigate the conservation of inflammatory routes. These proteins were divided into three groups: (i) cell marker proteins (colony-stimulating factor 1 receptor The gene sequences of the selected proteins were obtained from the gilthead seabream database using the whole-genome shotgun method (Pareek et al. 2011). The open reading frames of such genes were located using the ORFfinder software from the NCBI website (https:// www. ncbi. nlm. nih. gov/ orffi nder/), and an additional check was performed using NCBI Protein BLAST sequence alignment analysis (National Institutes of Health). The protein sequences obtained were used as input in the STRING database for the protein-protein interaction network functional enrichment analysis, comparing sequences between teleosts (superclass Actinopterygii). The percentages of identity of the most similar proteins and the e-values were used to investigate the conservation and interactions of gilthead seabream proteins. The protein names were derived automatically based on a cluster's consensus protein annotations taken from Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome (2021), Uni-Prot, Pfam, Simple Modular Architecture Research Tool (SMART), and InterPro. The predicted orthologous proteins were associated in networks, where nodes (proteins) and edges (protein-protein associations) denoted the number of predicted interactions Fish Physiol Biochem between them compared to interactions that would be expected for a random set of proteins of a similar size.

Design of the primers
The primers used were designed with the Thermo Fisher OligoPerfect™ tool according to the following criteria: (i) each oligonucleotide was composed of 20 nucleotides, (ii) the size of the amplicon had between 100 and 120 nucleotides, (iii) the guanidine-cytosine content was 55-60%, (iv) the semi-naturalization temperature (melting temperature) was as close to 60 °C as possible, and (v) primers that self-inhibit and form hairpins were avoided so as not to hinder the amplification reaction (Table 1).

Animals
Sixteen specimens (23.7 ± 7 g in weight and 11.7 ± 1 cm in length) of the seawater teleost gilthead seabream, obtained from a local farm (Mazarrón, Spain), were kept in recirculating seawater aquaria (250 L) in the Marine Fish Facilities of the University of Murcia (Spain) for a quarantine period of 1 month. The water temperature was maintained at 20 ± 2 °C with a flow rate of 900 L h −1 , 28% salinity, and a 12-h light and 12-h dark photoperiod. The water in the tanks was continuously aerated. The fish were fed a commercial diet (Perla MP, Skretting) in quantities amounting to 2% of their body weight daily and were subjected to 24-h fasting before the experiment. Commercial diet consisted of crude protein (48.5%), crude oils and fats (18%), crude ash (6.8%), crude fiber (2.2%), phosphorus (0.9%), calcium (0.9%), and sodium (0.3%). All experimental protocols were conformed to Directive 2010/63/EU and approved by the Ethics Committee of the University of Murcia.

Experimental design and sample collection
The fish were anesthetized with clove oil (20 mg L −1 ; Guinama), and subcutaneously injected in their left flank, in the middle part between the lateral line and the central zone of the anal fin. Two groups with two replicates (n = 4) were randomly established: (i) fish injected with 50 µl of phosphate-buffered saline (PBS; 11.9 mM phosphates, 137 mM NaCl, and 2.7 mM KCl; pH 7.4; Fisher Bioreagents) (control group) and (ii) fish injected with 50 µl of carrageenin in PBS (1%; Sigma-Aldrich) (carrageenin group). Twelve and 24 h after injection, two fish from each tank were randomly caught and sedated as previously described. Skin from the injected area was collected with a biopsy metal punch 4 mm in diameter (Stiefel). The skin samples were immediately stored in TRIzol Reagent (Invitrogen) at − 80 °C for gene expression analysis (Chomczynski 1993).

Gene expression analysis by real-time qPCR
Total RNA was extracted from 0.5 g samples of gilthead seabream skin using TRIzol, following the manufacturer's instructions, and quantification and purification were assessed using a Nanodrop® spectrophotometer (the 260:280 ratios were 1.8-2.0). Then, the RNA was treated with DNase I (Promega) to remove genomic DNA contamination and complementary DNA (cDNA) was synthesized from 1 µg of RNA using the reverse transcriptase enzyme SuperScriptIV (Life Technologies) with an oligo-dT 18 primer. The expression of the nominated genes (see Table 1) was analyzed by real-time qPCR with QuantStudio™ Real-Time PCR System Fast (Life Technologies). The reaction mixtures [containing 5 µl of SYBR Green supermix, 2.5 µl of primers (0.6 µM each), and 2.5 µl of cDNA template] were incubated for 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C, 1 min at 60 °C, and finally 15 s at 95 °C, 1 min at 60 °C, and 15 s at 95 °C. The gene expression was analyzed using the 2 −ΔCt method (Livak and Schmittgen 2001), which was performed as described elsewhere (Cordero et al. 2015). The specificity of the reactions was analyzed using samples without cDNA as negative controls. For each mRNA, gene expression was normalized with the geometric mean of ribosomal protein (s18), elongation factor 1-alfa (ef1α), and beta-actin (actb) RNA content in each sample. Gene names follow the accepted nomenclature for zebrafish (http:// zfin. org/). In all cases, each PCR was performed with triplicate samples.

Statistical analysis
The results were expressed as mean ± standard error of the mean (SEM). Data were analyzed by two-way ANOVA (followed by Tukey tests) to determine differences between experimental groups and each group with respect to time. The normality of the data Fish Physiol Biochem

Protein prediction analysis
To assess the conservation and functional interactions of the proteins involved in the inflammatory process, 40 gilthead seabream proteins were compared with protein sequences of other teleosts that appeared to match the input proteins in the STRING database ( Table 2). The large yellow croaker was selected for the analysis due to the higher number of proteins that matched the gilthead seabream protein sequence and the higher homology of its proteins (identity) with those of the gilthead seabream compared to the other teleosts. Gilthead seabream proteins were divided into three groups (cell markers and pro-inflammatory and anti-inflammatory proteins) according to their percentages of identity (Table 3). Cell marker proteins of the gilthead seabream sequence showed percentages of identity with that of yellow croaker ranging between 67.1 and 91.4%.
The identity of anti-inflammatory proteins was more homogeneous than that of pro-inflammatory proteins, with a range of identity between 59.9 and 91.3% (Table 3 family sequences showed a wider range, with NLRC3 (88.2%) and NLRX1 (86.9%) exhibiting higher identity percentages than NLRC5 isoforms 1 (70.3%) and 2 (59.9%).

Functional protein association network analysis
The functional association network analysis revealed that 34 nodes (functional proteins) selected from the gilthead seabream sequence interacted in large yellow croaker (Larimichthys crocea) ( Table 4). These 34 nodes represented 125 edges (specific and meaningful associations). The average node degree provided by STRING, which represents the number of interactions (at the score threshold) that a protein has on the average in the network, had a value of 7.35, while the clustering coefficient (a measure of how closely connected the nodes in the network are) had a value of 0.688 (Fig. 1).

Gene expression analysis
Real-time PCR was used to analyze the expression profile of 4 cell markers, 24 pro-inflammatory genes, and 12 anti-inflammatory genes in fish skin samples collected 12 and 24 h after carrageenin or PBS injection.

Cell markers
The expression of the csf1r gene (macrophage marker) 12 h after carrageenin injection was significantly downregulated in comparison to that of control specimens ( Fig. 2A). On the other hand, although no significant differences (p > 0.05) were observed either at 12 or at 24 h in the expression of phox22 and phox40 (acidophilic granulocyte markers) or mhciia (antigen-presenting cell marker) genes ( Fig. 2B-2D), an ascendant tendency in the expression of phox22 and phox40 was found in the carrageenin group 24 h after injection with respect to the control group.

Anti-inflammatory genes
The injection of carrageenin resulted in no significant changes in the expression of anti-inflammatory genes (il-10, tgf-β1, ctsd, ctsl, ctss, nlrc3, nlrc5 isoforms 1 and 2, nlrx1, ache, bche, and chrna7) with respect to the control group either at 12 or at 24 h (Figs. 7 and 8). All the results of the gene expression analysis are summarized in a schematic inflammatory model shown in Fig. 9.

Discussion
Although inflammation has been extensively studied in mammals, it is possible that not all the transduction pathways of this process have been conserved in the   (Byadgi et al. 2018;Savan and Sakai 2006). To our knowledge, this is the first study to investigate the conservation of inflammatory pathways in fish using the STRING database to identify functional associations of inflammationrelated proteins applicable to the gilthead seabream (Balasch and Tort 2019;Cordero et al. 2017;Sarropoulou et al. 2007;Silva-Marrero et al. 2017;Tsakogiannis et al. 2019). In our analysis, L. crocea, a species of particular interest for the marine aquaculture of the north-western Pacific, was the species whose protein sequence seemed to most closely match that of the gilthead seabream (Tan et al. 2019;Wang et al. 2016). It can be assumed that the analyzed proteins are highly conserved between the gilthead seabream and yellow croaker, since protein sequences with > 25% identity are considered orthologs (genes homologous by speciation) at the protein level and could be used to predict their molecular functions (Konaté et al. 2017). Furthermore, none of the studied sequences showed a percentage of identity less than 30%, the threshold below which protein structures degenerate (Gilson et al. 2017). Moreover, the 34 nodes obtained represented 125 edges, while the number of expected edges for a random set of proteins of the yellow croaker genome was 4. This high number suggests that the selected proteins are, at least partially, biologically connected (not necessarily physically bound but with some specific and meaningful associations), suggesting a shared function. These proteins are quite likely associated with inflammation (Postlethwait et al. 2000). It is noteworthy that the level of molecular conservation of cell markers CSF1R (also known as macrophage colony-stimulating factor receptor, MCSFR) and PHOX40, and PHOX22 between the gilthead seabream and large yellow croaker was significantly high (identity percentage close to 90%). These are important biomarkers for macrophages and acidophilic granulocytes (Grayfer et al. 2009;Meseguer et al. 1994). Interestingly, although the catalytic tyrosine kinase domains of CSF1R are highly conserved among vertebrates, it seems that csf1r genes exhibit a low sequence identity (especially in their extracellular domains), which could be reflected in the low conservation of the CSF1R ligand and the distinct macrophage functionality across divergent species (Grayfer et al. 2018). Likewise, the catalytic domains of PHOX22 and PHOX40, which are transmembrane and cytosol proteins, respectively, of the NADPH oxidase enzyme system of gilthead  (Belambri et al. 2018;Kawahara and Lambeth 2007;Meseguer et al. 1994;Sepulcre et al. 2002). In contrast, the basic pattern of MHC variation in mammals (Dijkstra et al. 2013;Yamaguchi and Dijkstra 2019) that has remained in fish with a high level of allelic/haplotype diversification explains the lower percentage of identity of MHC-II compared to the other cell markers in our study. Regarding pro-inflammatory proteins, we found a high conservation level of proteins related to the NF-κB family, consisting of five polypeptide subunits-RelA (p65), RelB, C-Rel, NF-κB1 (p50/ p105), and NF-κB2 (p52/p100). This could be due to the key role that this transcription factor plays in the activation and release of pro-inflammatory cytokines, as well as its implication in the expression of pro-survival molecules involved in the regulation of important processes, such as the immune response and homeostasis maintenance (Etemadi et al. 2015;Gugasyan et al. 2004;Mulero et al. 2019;Napetschnig and Wu 2013). This is also supported by the relatively highly conserved adaptors and regulatory proteins (STAT3, TRAF6, IκBKG, MYD88, and IRAK1), which might be critical for maintaining a minimal threshold of NF-κB signalling and whose activity could also be regulated by intracellular homeostatic processes, such as oxidative stress, metabolite changes, and DNA damage (Hinz et al. 2010;Linares et al. 2013;Matsuzawa et al. 2005;Moscat and Diaz-Meco 2016;Rezaeian et al. 2017;Tang et al. 2013;Tzeng et al. 2013;Zhang et al. 2016).
Likewise, the high percentage of identity (around 85%) of anti-inflammatory proteins, such as cholinergic proteins and cathepsins, could also be related to their regulatory function. For instance, facing an inflammatory response, the nervous system of higher organisms can release acetylcholine through the efferent vagus nerve, which is activated by AChE and BChE, and inhibit the production of pro-inflammatory cytokines in macrophages via activation of the α7-nicotinic homopentameric receptor (α7nAChR) (Zila et al. 2017). It has been observed in human cell lines that NOD-like receptors, such as NLRC3, may also downregulate the expression of NF-κB transcription factor by TRAF6 ubiquitination, while mitochondrial NLRX1 and NLRC5 limit the NF-κB activation by interfering with the interactions of the TRAF6-IKK (Iκ-B kinase) protein complex and the IKKα protein, respectively (Cui et al. 2010;Schneider et al. 2012;Xia et al. 2011). Cathepsins, for their part, not only are lysosomal endonucleases that negatively modulate inflammation but are also involved in cell survival and apoptosis, participating, among other functions, in neutrophil recruitment (CTSL), caspase activation (CTSD), and process and maturation of the associated invariant chain of MHC-II (CTSS) Li et al. 2015).
It is important to emphasize the ancient genome duplication that occurred in teleosts, whereby an estimated 30% of genes were duplicated. This could explain the main protein sequence differences found in our assay (Postlethwait et al. 2000). It could also support the high conservation of molecules involved in the regulatory pathways of the immune response (such as the NF-κB family transcription factors). On the other hand, secondary mediators (such as proinflammatory and anti-inflammatory cytokines and proteins of the TLR and NOD families and the TNR superfamily) could vary structurally or sequentially depending on the organism, preserving their functions to a greater or lesser extent in most teleost fishes as paralogues (multiple copies) (Mulero et al. 2019;Secombes et al. 2011;Wang and Secombes 2013). This also explains the presence of two isoforms of NLRC5 in the gilthead seabream, whose protein sequence seems to coincide to some degree with the NLRC5 protein sequence of L. crocea, suggesting that NLRC5 isoform 1 of the gilthead seabream could be more conserved than isoform 2.
To validate the in silico results obtained from the STRING database and to confirm the conservation of the molecules involved in the inflammatory response of the gilthead seabream, an in vivo analysis of the gene expression was performed using carrageenin as a stimulus of inflammation according to the schema shown in Fig. 9. The concentration of injected carrageenin (50 µl of 1% carrageenin) and the sampling times (12 and 24 h after injection) were carefully chosen considering the available data from studies on mammals (Levy 1969;Morris 2003;Winter et al. 1962). Nonetheless, although proteins seem to be relatively highly conserved, some key considerations are important for understanding the functional differences between studied organisms and the effects of carrageenin on the gilthead seabream. For instance, unlike in mammals, leucocytes in fish should migrate from the head kidney (the main hematopoietic organ) to the injured site when an insult is perceived by local cells (Meseguer et al. 1995), such as the presence of carrageenin in this study. Therefore, infiltrating acidophilic granulocytes of the gilthead seabream, which are functionally equivalent to mammalian neutrophils (Meseguer et al. 1994;Sepulcre et al. 2002), should be the first cells to be recruited to the inflammation site, followed by macrophages (Kolaczkowska and Kubes 2013;Medzhitov 2008;Nathan 2006;Nguyen-Chi et al. 2015). The observed downregulation of csf1r expression 12 h after carrageenin injection suggests a Fig.1 Representation of the node network (selected proteins) of Larimichthys crocea generated from the Sparus aurata protein sequence using the STRING database (see protein annotation in Table 3). The different edge colors indicate the type of evidence for each interaction: cyan = curated STRING data-base; purple = experimentally determined; green = gene neighborhood; blue = gene co-occurrence; light green = text mining; black = co-expression; lavender = protein homology (the types of evidence are not mutually exclusive) Fish Physiol Biochem  Fig. 2 Relative expression of A csfr1, B phox22, C phox40, and D mhciia in gilthead seabream injected with PBS (control, white bars) and carrageenin (1%, black bars) analyzed in skin samples collected 12 and 24 h after injection. The bars represent the mean ± standard error of the mean (n = 4). The asterisks denote significant differences between the carrageenin and control groups (two-way ANOVA; p < 0.05)   (Kaneda et al. 1991;Mitsuyama et al. 1982;Phelps and Neely 2007), which lends weight to the hypothesis of a cytotoxic effect of carrageenin on gilthead seabream macrophages. Furthermore, the observed increase tendency of phox22 and phox40 at both 12 and 24 h after carrageenin injection suggests the activation or increased recruitment of acidophilic granulocytes to the inflammation site at that time. However, further studies are needed to elucidate the effect of carrageenin at the subcellular level and confirm these hypotheses (Belambri et al. 2018). Otherwise, and regarding the subunits of NF-κB transcription factor in mammals, it would be important to consider that RelA:NF-κB1 heterodimer is the principal product that induces the canonical signalling pathway of inflammation. RelB:NF-κB2 constitutes the non-canonical or alternative activation pathway, leaving c-REL subunit as a secondary molecule in this process (Mulero et al. 2019). Interestingly, the upregulated expression of c-rel at 24 h found in our study after carrageenin injection could suggest an alternative pathway of activating inflammation in fish. Thus, c-REL homodimers could combinate activating the transcription of inflammation-related genes, since c-REL heterodimers with other NF-κB components have not been documented (Mulero et al. 2019). In addition, the trigger of a faster inflammatory response than in mammals could also explain that the expression of cytokine genes was not significantly altered in the time tested. Taken together, these results suggest that the activation of these genes might be the first step in the termination of the innate immune response in order to restore homeostasis (Gugasyan et al. 2004).

Conclusion
These results indicate that the molecules involved in the inflammatory process, especially the proteins involved in the regulation of the inflammatory Relative expression of A nlrc3, B isoform 1 of nlrc5, C isoform 2 of nlrc5, D nlrx1, E ache, F bche, and G chrna7 in gilthead seabream injected with PBS (control, white bars) and carrageenin (1%, black bars) analyzed in skin samples collected 12 and 24 h after injection. The bars represent the mean ± standard error of the mean (n = 4). No significant differences were obtained neither between the carrageenin and control groups nor between the sampling time points (two-way ANOVA; p < 0.05) Fish Physiol Biochem response, are mostly conserved among fish species, but their main functions can vary depending on the species due to the different orthologues present in each. Furthermore, the inflammation induced by subcutaneous carrageenin injection (at least at the tested dose) in the gilthead seabream abates within hours. The mechanisms necessary to terminate the inflammation response and recover skin homeostasis are activated between 12 and 24 h after injection. The in silico and gene expression analyses performed in this study might contribute to the identification of the main mechanisms of the acute inflammatory response. Our results have implications not only for applications in the aquaculture sector but also for basic research. Future studies using this characteristic pro-inflammatory mucopolysaccharide are needed to elucidate the complex mechanism of inflammation related to diseases or dietary supplementation in fish.
Author contribution MAE conceptualized the study and MAE and FAG reviewed the research proposal. JCCS and JML conducted the experimentation and sampling. JCCS and JML processed the samples. JCCS and FAG analyzed the data. JCCS and JML wrote the manuscript and MAE and FAG reviewed it. MAE got the financing. All the authors participated in giving feedback on the manuscript. All authors have read and approved the final manuscript.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. This work was supported by the Spanish Ministry of Economy and Competitiveness (Grant no. AGL-2017-83370-C3-1-R) co-funded with Fondos Europeos de Desarrollo Regional (FEDER) and Fundación Séneca de la Región de Murcia (Grupo de Excelencia 19883/GERM/15).

Availability of data and material
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Code availability Not applicable.

Declarations
Ethics approval/declarations The experiment was conducted at the Marine Fish Facilities at the University of Murcia (Spain). All experimental protocols were approved by the Ethical Committee of the University of Murcia (permit number CEEA 357/2017) following the guidelines of European Union for animal handling (2010/63/EU).

Consent to participate
All authors consent the participation.

Conflict of interest
The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.